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Single-shot measurement of the Josephson charge qubit 
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We demonstrate single-shot readout of quantum states of the Josephson charge qubit. The quan- 
tum bits are transformed into and stored as classical bits (charge quanta) in a dynamic memory 
cell - a superconducting island. The transformation of state |1) (differing form state |0) by an extra 
Cooper pair) is a result of a controllable quasiparticle tunneling to the island. The charge is then 
detected by a conventional single-electron transistor, electrostatically decoupled from the qubit. We 
study relaxation dynamics in the system and obtain the readout efficiency of 87% and 93% for |1) 
and |0) states, respectively. 

PACS numbers: 03.67.-a, 74.50.+r, 85.25.Cp 
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It has been recently realized that Josephson junctions 
can be used for building quantum bits (qubits) and inte- 
grated quantum computer circuits controlled by external 
electrical signals J^i^. After the first experiments on sin- 
gle Josephson qubits 0, S S IE 1 an important step to- 
wards the integration has been made: coherent control of 
two qubits and conditional gate operation have been ex- 
perimentally demonstrated for two electrostatically cou- 
pled charge qubits |2i,||9||- However, in these experiments, 
individual probabilities of each qubit averaged over all 
states of the other qubit were measured fl^. To directly 
measure multi-qubit states, one must be able to readout 
each qubit after every single-shot coherent state manipu- 
lation. The single-shot readout is of great importance, for 
instance, for quantum state tomography, quantum state 
teleportation, quantum cryptography [Tl|. Without the 
single-shot readout, algorithms that give non-unique so- 
lutions can not be utilized. 

To readout single quantum states of the Josephson 
qubits (in particular, flux qubits) through the phase de- 
gree of freedom, a few circuits, measuring switching event 
from the supercurrent state to the finite-voltage state 
were implemented 0, IE IE 0- charge type of qubits, 
it is straightforward to measure a charge quantum in- 
stead of the flux quantum fl3| . For the single-shot charge 
readout, a radio- frequency single-electron transistor [l3l | 
electrostatically coupled to the qubit was proposed as a 
detector of the charge stat es Il 4[ll5ll . Although this ap- 
proach works in principle |l6l Il7l| , the single-shot read- 
out has not yet been realized. In this work, we demon- 
strate an operation and study mechanism of novel read- 
out scheme that allows to perform highly efficient single- 
shot measurements, with suppressed back-action of the 
measurement circuit on the qubit. 

A scanning-electron micrograph of our circuit is shown 
in Fig.nja) The device consists of a charge qubit and 
a readout circuit. The qubit is a Cooper-pair box (with 
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its effective capacitance to the ground Cb ~ 600 aF) cou- 
pled to a reservoir through a Josephson junction with 
the Josephson energy Ej k, 20/zeV. The reservoir is a 
big island with about 0.1 nF capacitance to the ground 
plane and galvanically isolated from the external envi- 
ronment. The qubit states are coherently controlled by a 
non-adiabatic control pulse, yielding a superposed state 
of |0) and |1). The readout part includes an electrome- 
ter which is a conventional low-frequency single-electron 
transistor (SET) {Cs « 1000 aF) and a charge trap (Ct « 
1000 aF) placed between the qubit and the SET. The 
trap is connected to the box through a highly resistive 
tunnel junction {Rt w 100 MD) and coupled to the SET 
with a capacitance Cst ~ 100 aF. The use of the trap 
enables us to separate in time the coherent state manip- 
ulation and readout processes and, in addition, the qubit 
becomes electrostatically decoupled from the SET. The 
qubit relaxation rate induced by the SET voltage noise 
is suppressed by a factor of [CbtCst/CtCbY ~ 3 x 10~^, 
where Cbt « 30 aF |3. The coupling strength can be 
made even weaker, if the unwanted box-to-trap capaci- 
tance Cbt is further decreased. 

The operation of the circuit can be described in the fol- 
lowing way. During the qubit manipulation, the trap is 
kept unbiased prohibiting charge relaxation to the trap. 
Once the control pulse is terminated, the readout pulse 
(see the inset of Fig. Q^a)) is applied to the trap. The 
length and the amplitude of the readout pulse are ad- 
justed so that if there is an extra Cooper pair in the box 
after termination of the control pulse, it escapes to the 
trap through a quasiparticle tunneling with a high proba- 
bility. After the charge is trapped, it remains in the trap 
for a long time (a reverse trap-to-box charge relaxation is 
suppressed due to the superconducting energy gap 2A) 
and is measured by a low-frequency SET. 

The Hamiltonian of the two-level system of the qubit 
in the charge basis |0) and |1) (without and with an ex- 
tra Cooper pair) is = i7b(0, Qb)|0)(0| + C/fc(2, Qh)|l)(l| - 
1/2 Ej (|0)(1| -I- |1)(0|) (we define an electrostatic energy 
of island k as Uk{Nk,Qk) = {NkC - Qkf /2Ck, where k 
is either 6 or i indicating box or trap island, respectively, 
Nk is an excess electron number, and Qk is a gate in- 
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FIG. 1: (a) Scanning electron micrograph of the device. The 
aluminum structure is deposited on top of a thin Si3N4 in- 
sulating layer (0.4 ^m) above a gold ground plane. The de- 
vice consists of a Cooper pair box, a reservoir, a trap and 
a measurement SET. The dc ("box" and "trap") and pulse 
(control and readout) gates control potentials of the islands. 
Pulse operation is schematically represented in the inset, (b) 
Stability diagram of the SET coupled to the trap. Open cir- 
cles mark positions of the SET quasiparticle current peaks on 
Vgs — Vgt plane (dc gate voltages of the SET and the trap). 
Pairs of numbers {Nt,Ns) designate the trap - SET ground 
state charge configuration in each cell bounded by the SET 
peaks and solid lines. Dashed, dashed-dotted and dotted lines 
indicate positions of the SET peaks for 0, 1 or 2 additional 
electrons in the trap, respectively, (c) A typical time-trace 
of the SET current (lower panel) together with the readout 
pulse sequence (upper panel). Negative switches on the lower 
curve correspond to the detected charge of the trap. Digits 1 
and mark readout bits. 

duced charge in the island. Starting at QboiQbo < Qbi), 
where E » Ej {AE = Ut{2, Qt) - Ut{0, Qb)) we let the 
system relax to the ground state, which is nearly pure 
charge state |0). Then we instantly change the eigenba- 
sis for a time tc by applying a rectangular control pulse, 
which brings the system to Qbi- If Qbi is a degeneracy 
point (Ai? = 0), the final state of the control pulse ma- 
nipulation is |0) cos + |l)sinwjtc/2 {ojj = Ej/h), 
therefore, after the pulse termination, the state |1) is re- 
alized with a probability of cos^ ujjtc/2. 

Fig. n^b) shows an experimentally measured stability 
diagram: SET current peak positions as a function of 
trap and box gate voltages of the SET and the trap. 
By setting the box and trap gates to one of the points 
Nto, Nti or Nt2, we can detect if the trap has 0, 1 or 
2 additional electrons. In our measurements, the SET is 
usually adjusted to Nto- To readout the qubit, the trap is 
biased by the readout pulse of typical length U ~ 300 ns 
and amplitude AQt — 3.5e (AQf = Qti — Qto), applied to 
the readout gate, letting an extra Cooper pair of the state 
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FIG. 2; (a) Coherent oscillations measured by averaging 
over many events of the single-shot measurements, (b) P 
versus tc measured at AQj, = 0.84e (= QbA), where visi- 
bility is the highest, (c) P versus tc measured at Qb = 0.75e 
(= Qbs), where the oscillations are the longest lasting (degen- 
eracy point). Dashed envelops correspond to the exponential 
decay with the decay time of 5.8 ns. 



|1) escape to the trap through a quasiparticle tunneling 
and switching off the SET current (the SET peak position 
is shifted to the position of the dashed line). 

The curve in the upper panel of Fig. ^c) indicates a 
readout pulse sequence. The curve on the lower panel 
demonstrates a typical time-trace of the SET current. 
Negative switches on the curve of the lower panel com- 
ing synchronously with the readout pulses are counted 
as charge detection events. For the studied device, the 
lifetime of the trapped charge is typically about 300 /iS, 
therefore, normally used repetition time of 2 ms is suf- 
ficient for practically complete trap resetting. We count 
the number of detected switches to, with the total num- 
ber of shots Tit. 

An experimentally obtained probability of the charge 
detection P = m/rit (normally, rit = 327 is used per 
one experimental data point) as a function of the con- 
trol pulse length t,. and the amplitude AQf, {AQt = 
Qti — Qto) is shown as a two-dimensional plot in Fig.j^Ja). 
We define the pulse with AQb = 0.84e (= QbA) and 
tc — 120 ps, when P reaches maximum, as a 7r-pulse. 
Fig. Elb) demonstrates coherent oscillations as a func- 
tion of ic measured at AQb a- As shown by the vertical 
arrowed line, the visibility here reaches 0.64, while the 
longest lasting oscillations shown in Fig. [21[c) are found 
at AQb = 0.75e (= Qbs)-, (the phase decoherence is ex- 
pected to be the weakest at the degeneracy point). We 
beHeve that QbA QbB due to the control pulse dis- 
tortion because of limited frequency bandwidths of the 
transmission lines and the pulse generator. 

Curves with open and closed dots in Figs. |3fa-b) rep- 
resent P as a function of the readout pulse amplitude 
AQt with and without control 7r-pulses, respectively. The 
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FIG. 3: P versus readout pulse amplitude AQt with and 
without TT-puIse (open and closed cycles, respectively), (a) 
P measured at Nto (see Fig. ^b)). (b) P measured at Nt2- 
One may distinguish three different regions on the plots: (I) a 
finite probability for the excited state detection is only at Nto 
position; (II) finite probabilities for the excited state detec- 
tion are at A'^to and Nt2 positions; (III) switches are detected 
even if the qubit is in the state |0). Threshold amplitudes 
QtA, QtB, Qtc are derived for the box-to-trap relaxation pro- 
cesses (2,0) ^ (1,1), (1,1) ^ (0,2) and (0,0) ^ (-1,1), 
respectively. 



probabilities measured at Nto are shown in Fig. Ola), 
while the probabilities measured at Nt2 (where positive 
switches are counted) are shown in Fig. |2tb). One may 
divide the plots into three regions marked by I, II and 
III differing from each other by the counting character- 
istics. Based on the data from these plots, we suppose 
that when the qubit is in the state |1), one quasiparticle 
tunnels from the box to the trap in the process (2, 0) 
(1, 1) within the region I {{Nb,Nt) represents the box- 
trap quasiparticle configuration); quasiparticles tunnel to 
the trap in the process (2, 0) (1, 1) — > (0, 2) within 
the region II; a quasiparticle tunneling process becomes 
possible even for the state |0) in the region III. 

The quasiparticle tunneling is energetically feasible in 
the process {Nb,Nt) {Nt - l,Nt + I), when the fol- 
lowing condition is satisfied Ub{Ni,, Qto) + Ut{Nt,Qti) > 
Ub{Nb - 1, Qbo) + Ut{Nt + 1, Qti) + 2A 19] (we neglect 
the interaction energy term which is as small as the inter- 
island coupling strength). Substituting an explicit ex- 
pression for the energies, one may find the necessary 
trap readout pulse amplitude for the quasiparticle es- 
cape AQf > (3j + A^fe + ?7(2-iVb)e, where 77 = Ct/Cb and 
Q[ = e/2-77(3e/2-gbo) + 2ACt/e-Qto- We define three 
threshold amplitudes of AQt, at which the following pro- 
cesses become possible: (2, 0) (1, 1) at AQt > Qf^, 
where QtA = QJ ; (1, 1) ^ (0, 2) at AQt > QtB-, where 
AQts = Q't + (1 + ri)e- (0, 0) ^ (-1, 1) at AQt > Qtc, 
where Qtc = Q't + 2776. The threshold amplitudes cal- 
culated with QtA = 0.6e (taken from the first step of P 
curve on Fig. O^a)) and i] = 1.67 are QtB = 3.3e and 
Qtc — 3.9e and shown on the Figs.|3ta-b) by the dashed 
vertical lines. Note that the probability for the trap to 
have a charge at AQt > 4.5e reaches 1 in Fig. ISJa). This 
proves that once the charge is trapped it is detected with 
100% probability. For the highest efficiency, we set the 
operation point of the SET to A'to, where either one or 
two trapped quasiparticles yield a negative switch of the 



FIG. 4: (a) P as a function of delay between the control tt- 
pulse and the readout pulse td- A solid curve is a fitting by an 
exponent with a decay rate fi = (220 ns)~^. (b) A probability 
of |l)-states detection created by the -pulse as a function of 
the readout pulse length tr . A solid curve is a result of fitting 
m(t) using Eq. Q normalized by nt with fitting parameters 
no and a. {n*)/nt = 0.87 and a = (37 ns)~^). The inset shows 
probability without 7r-pulses. P is fitted by 1 - exp(— 7tr) 
with 7 = (4.1 ^s)~^. 

SET current and the pulse amplitude to AQt w 3.5e, at 
which two quasiparticles may escape to the trap. 

Figs, ^a-b) demonstrate time relaxation dynamics of 
the qubit states. Fig. 01a) shows a probability P to find 
an extra charge in the trap, when time delay td is in- 
troduced between the control 7r-pulse and the readout 
pulses. The exponential decay of P may be explained by 
tunneling to the reservoir (presumably via energetically 
feasible Cooper pair tunneling (2, 0) — » (0, 0)) because 
alternative quasiparticle relaxation to the trap through 
the high resistive junction is blocked by 2A when the 
trap is not biased. The relaxation rate to the reservoir 
found from the fitting (sohd curve) is /3 = (220 ns)~-^. 
Fig- Etb) shows relaxation dynamics of the state |1) as 
a function of the readout pulse length tr {td ~ 0). This 
relaxation is mainly determined by quasiparticle decay 
to the trap with a rate a (a 3> /?). Additionally, the 
inset of Fig.^fb) shows dynamics of |0)-state relaxation 
("dark" switches) to the trap. These "dark" switches can 
be presumably described by the process (0, 0) (-2, 2), 
with a weak relaxation rate /i — (4100 ns)~^ derived from 
fitting the data by 1 - exp(— 7/:^) (solid curve). 

Let us consider the relaxation dynamics in more de- 
tails. The number of excited states, n*, decreases within 
the time interval [t, i -f dt] as dn*{t) — —an*{t)dt — 
(3n*{t)dt [231 . The number of states in (0, 0) configu- 
ration, n(t), changes, in turn, as dn(t) = ~jn(t)dt + 
j3n*{t)dt. We may also write an expression for the num- 
ber of events in which the trap is found to be charged: 
dm(t) = n*(t)dt + n{t)dt. Solving these equations with 
the initial conditions n*{Q) — n^, ^(0) — nt — Hq and 
to(0) =0 we find 

mit) = nt[l - + *:i2l[e-7t _ e-("+«*]. (1) 

a + p — "f 

We fit the data of Fig. 01b) by a curve of P = m{t)/nt 
with m{t) taken from Eq. |^ with two fitting parame- 
ters a and Uq. The fitting gives a = (37 ns)~^ |23 and 
n^/rit = 0.84, implying that the efficiency of |0)-to-|l) 
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conversion by the control 7r-pulse is 84%. 

If our readout pulse length t ^ tr satisfies the condi- 
tions 7 ^ <^ a + P, then Eq. ^ can be simpHfied 
to 



m{t) 



-f3' 



(2) 



Using Eq. JSJ, one may estimate an efficiency of the 
single-shot readout. We introduce a probability Py{x) 
of finding the trap charged (j/ = 1) or uncharged (y = 
0), when the qubit is in |x)-state {x is either or 1). 
According to the definition of Py{x), Po{0) + Pi(0) = 1 
and Po(l) + -Pl(I) = 1- The total number of detected 
events expressed in terms of these probabilities may be 
written as m = nQPi(l) + (rif — nQ)Pi(O). Comparing the 
latter expression with Eq. Q we find 



a + j3jtr 
a + P 

itr- 



(3a) 
(3b) 



Confirming that our readout pulse length tr ~ 300 ns 
fulfills the necessary condition for Eqs. Q, 7 ^ ^ 
q; -I- /?, we directly find from Eqs. Q that the probabil- 
ity of detection of the state |1) is Pi(l) = 0.87 and the 
state |0) is Po(0) = 0.93 (Pi(0) = 0.07). The readout ef- 
ficiency can be further improved by optimizing the relax- 
ation rates. The derived probabilities are consistent with 



the mean probability of the oscillations at the degener- 
acy point, (P> - [(nS)Pi(l) + (nt-(n5))Fi(0)]/nt = 0.47 
(horizontal line in Fig. QJc)), where (no/rit) = 0.5. 

Finally, we would like to briefly discuss new possibili- 
ties and potential applications, which this scheme opens. 
The trap island works as a dynamic memory cell, which 
stores the multi-qubit state in a classical bit. Therefore, 
the data can be further processed in a classical way. For 
example, implementing sequential readout of many cells 
connected in series, similarly to that of charge-coupled 
devices will allow to reduce a number of output circuit 
channels. 

One of straightforward applications of the described 
device, staying outside quantum computation direction, 
is a tunable single-photon detector in a centimeter wave- 
length range. A weak microwave radiation may reso- 
nantly excite the box near the degeneracy point from 
the ground to the excited state with a finite probabil- 
ity. The excited state may then be converted into the 
charge state |1) by an adiabatic sweep of the control gate 
voltages (a relatively slow gate voltage change forces the 
system to move along eigenenergy bands) and detected 
by the single-shot measurement circuit. 

As for the nearest plans, it is straightforward to use the 
single-shot readout to for measuring double qubit system 
|@i [9| . This will allow us to demonstrate entangled states 
and controUed-NOT gate operation. 

We thank S. Lloyd for valuable discussion. 
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